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摘要 :分 别 以 18SrDNA 的 V4 区 和 V9 区 为 目标 基因 ,采用 高 通 量 测序 平台 和 生物 信息 学 方法 ,分 析 海 水 样品 中 微型 和 微微 型 浮 
游 植物 多 样 性 。 利 用 在 线 分 析 软 件 对 VACE/R) , V9( F/R) flü CA(E/R)3 对 引物 的 敏感 性 人 特异 性 进行 了 评估 和 比较 ,发 现 自行 
设计 的 引物 V4(F/R) 对 真 核 池 类 的 扩 增 特异 性 高 于 V9( F/R) 和 C4( F/R)。 席 通 量 测序 结果 显 东 ,检测 样品 平均 获得 68834 
条 原始 序列 ,高 质量 数据 占 99% 以 上 ,获得 基因 注释 的 序列 数 达 94% 以 上 。3 对 引物 V4(FAR) 、V9(F/R)、C4(F/R) 鉴定 的 平 
均 微型 /微微 型 浮游 植物 OTUs 数 分 别 为 78 .42 .58 ,引物 V4(F/R) 鉴 定 效率 由 对 较 高 ,同时 对 细小 微 胞 藻 ( Micromonas pusilla) 、 
(A ^F- fs ER HE. Ostreococcus tauri) , $8 ER $E ( Pycnococcus provasolii ) WE 4 EREE Aureococcus anophagefferens ) , 2j W 5t 25 W 
C Heterosigma. akashiwo ) 等 优势 种 检 出 频率 高 于 引物 V9( F/R)。 相 对 妃 发 表 的 2 对 引物 ,设计 的 引物 V4( F/R) 对 海洋 微型 和 微 
微型 藻类 多 样 性 检测 更 为 高 效 。 
关键 词 : 微 型 浮游 植物 ; 微微 型 浮游 植物 ; 高 通 量 测序 ; 18S rDNA X V4; 可 变 区 V9 


si 


Optimization of high-throughput sequencing primers for nanophytoplankton and 


picophytoplankton in environmental samples 
LIU Weidong"^" , SONG Lun'^, WU Jing ^ 


1 Liaoning Ocean and Fisheries Science Research, Institute , Dalian 116023, China 


2 Liaoning Ocean Environment Monitoring Station , Dalian 116023, China 


Abstract: In this study, nanophytoplankton and picophytoplankton diversity in seawater samples was analyzed using a high- 
throughput sequencing platform and a series of bioinformaties tools, based on the V4 and V9 region of 18S rDNA as the 
target gene. High-throughput sequencing, which is considered as one of the most important tools in genomics research, is 
widely applied in the field of marine nanophytoplankton and picophytoplankton diversity studies. We successfully obtained a 
pair.of nanophytoplànkton and picophytoplankton PCR primers V4( F/R) by analyzing the nucleic acid database and using a 
series of bioinformatics tools. Two pairs of universal primers were also selected for comparative analysis, which amplified 
variable region V4 and V9 of the small subunit nuclear ribosomal DNA ( SSU nrDNA ) .The sensitivity and specificity of PCR 
primers V4( F/R), V9( F/R), and CA( F/R) were also evaluated and compared using online bioinformatics software. The 
results showed that the amplification specificity of primer pair V4( F/R) was better than that of V9( F/R) and CA( F/R) in 
eukaryotic algae. High-throughput sequencing results showed that 68834 raw tags were amplified by the primers, 9996 of 
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which were effective tags. Sequences of more than 94% of the effective tags were identified by Ribosomal Database 
ProjectClassifier, among which 308 operational taxonomic units ( OTUs) of one sample were used for further analysis. The 
average numbers of nanophytoplankton and picophytoplankton OTUs amplified by V4( F/R), V9(F/R), and C4( F/R), 
were 78, 42, 58 ,respectively. The primer pair V4( F/R) was found to have higher sensitivity and specificity for amplifying 
nanophytoplankton and picophytoplankton, including Micromonas pusilla, Ostreococcus tauri, Pycnococcus provasolii, 
Aureococcus anophagefferens , and Heterosigma akashiwo. The V4 region from the environmental eukaryotic 188 rDNA gene 
could be suitable for high-throughput sequencing technology , and it was also a good target gene formarine nanophytoplankton 
and picophytoplankton identification. This study demonstrates the use of a simple, rapid, high sensitivity, and' low-cost 
technology to explore marine nanophytoplankton and picophytoplankton diversity. Moreover, it also provides a reference for 


the early warning and control of brown tide disasters. 


Key Words: nanophytoplankton; picophytoplankton; high-throughput sequencing; variable region V4. of 188 rDNA; 


variable region V9 


海洋 微型 和 微微 型 浮游 真 核 藻类 具有 生长 快 、 密 度 高 分布 广 等 特点 ,是 海洋 初级 生产 力 的 重要 贡献 者 ， 
同时 也 是 赤潮 , 褐 潮 暴 发 的 致 灾 种 ' 中 。 由 于 微型 和 微微 型 藻类 个 体 微 小 ,缺乏 形态 学分 类 指标 、 多 数 难 以 分 
离 培养 .无 法 大 量 单纯 富 集 ,所 以 传统 的 分 类 学 方法 鉴定 较为 困难 。 分 对 生物 学 的 快速 发 展 ,为 微型 和 微微 型 
藻类 的 多 样 性 研究 带 来 了 新 突破 :4 。Handelsman 等 5 于 1998 年 正式 提出 了 宏基 因 组 学 ,主要 以 环境 样品 
中 的 微生物 群 基因 组 为 研究 对 象 ;以 功能 基因 筛选 测序 分 析 为 研究 手段 ;以 微生物 多 样 性 种群 结构 、 系 统 进 
化 及 与 环境 之 间 的 关系 为 研究 目的 。 其 研究 对 象 已 从 最 初 的 土壤 微生物 发 展 到 水 体 微 生物 和 海 样 浮游 生 
物 !5 。 宏 基因 学 主要 的 研究 手段 之 一 是 根据 核糖 体 基 因 AxDNA 数据 库 进 行 引物 设计 ,通过 系统 发 生 学 分 析 得 
到 该 环境 中 生物 的 遗传 多 样 性 和 分 子 生态 学 信息 ;1 。 核 糖 体 基因 rDNA 序列 进化 变异 频率 缓慢 ,是 生物 群落 
系统 进化 分 类 研究 主要 基因 之 一 , 它 包括 糖 体 小 吏 基 基因 18S rDNA 、5.8S rDNA .大 亚 基 基因 28S rDNA 和 转 
录 间 隔 区 ITS。 其 中 核糖 体 小 亚 基 基因 18S,rDNA 罕 在 于 所 有 真 核 生物 中 ,是 同时 具有 信息 和 功能 两 种 作用 
的 核糖 体 编码 基因 ' ,其 数据 库 丰 富 基因 序列 种 类 广泛 ,对 引物 设计 构建 物种 亲缘 关系 等 方面 都 具有 不 可 
替代 的 优势 。18S rDNA 由 保守 区 和 可 赛区 交 蔡 排列 组 成 ,具有 9 个 可 变 区 '” ,可 作为 种 属 鉴定 的 分 子 标记 。 

Moon-van der Staay 和 LópezGarcía 通过 构建 18S rDNA 克隆 文库 和 分 子 系统 树 ,分 别 研 究 了 太平 洋 近 赤道 
海域 表层 海水 和 南极 深海 中 微微 型 真 核 浮游 生物 的 多 样 性 2 。 在 我 国 邵 鹏 、 袁 杰 等 人 分 别 对 厦门 西海 域 
和 南沙 群岛 海域 微型 和 微微 型 序 游 藻类 建立 了 18S rDNA 克隆 文库 ,分 析 了 该 水 域 浮 游 藻 类 的 多 样 性 站。 
自 2005 年 起 ,高 通 量 测序 技术 发 展 飞 速 ,由 于 该 技术 具有 简单 快速 .测序 通 量 高 .错误 率 低 和 成 本 低廉 等 特 
点 ,加 速 了 宏基 因 组 学 的 发 展 ] 。Stoeck 和 Behnke 等 人 对 18s rDNA 的 可 变 区 在 真 核 生物 多 样 性 研究 中 的 作 
用 进行 了 分 析 8 汐 基因 标记 的 选择 提供 了 参考 5 7 V4 区 和 vo 区 是 目前 广泛 应 用 于 真 核 生物 多 样 性 研究 
的 目标 基因 上 Amaral-Zettler 通过 以 18s rDNA 可 变 区 V9 为 目标 基因 设计 引物 ,利用 高 通 量 测序 技术 研究 了 海 
TETTE EI Z PEE! 。Stoeck 等 对 以 18s rDNA 的 V4 和 V9 区 为 目标 基因 ,结合 高 通 量 技术 获得 了 挪威 峡 湾 
TC Ug ^E DURS Ae RETE 。 

由 于 不 同 引 物 对 环境 中 生物 多 样 性 的 检测 效率 差异 较 大 ,针对 微型 和 微微 型 浮游 植物 ,本 研究 分 别 以 
18SrDNA 的 V4 区 和 V9 区 为 目标 基因 设计 和 选择 引物 进行 扩 增 ,对 引物 的 敏感 性 和 特异 性 进行 评估 和 比较 
分 析 ,采用 高 通 量 测序 平台 和 生物 信息 学 方法 ,应 用 于 辽东 湾 海域 进行 微型 和 微微 型 浮游 植物 检测 验证 ,筛选 
高 效 检测 引物 。 


1 材料 和 方法 


1.1 引物 设计 
本 文 针 对 目前 已 知 的 微型 /微微 型 真 核 藻类 以 及 黄 渤海 近 岸 赤潮 生物 种 类 设计 了 两 对 引物 。 首 先 于 
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NCBI 数据 库 调 取 相关 序列 ,其 中 包括 微微 型 真 核 藻 类 45 种 中国 近 岸 赤潮 生物 71 种 及 常见 真 核 藻类 73 种 ， 
涵盖 8 个 门 、18 个 纲 、179 条 序列 。 利 用 软件 MEGA4 对 序列 进行 多 重 序列 比 对 ,找到 对 应 的 V4 区 ,同时 结合 
Hiseq2500 PE250 测序 技术 要 求 进 行 引物 设计 ,片段 长 度 为 300 一 450bp ,为 增加 测序 结果 的 准确 性 ,充分 利用 
测序 技术 ,以 V4 区 为 核心 区 及 两 侧 的 保守 区 设计 引物 ,获得 两 对 引物 V4( F/R) ,使 用 FPCR 软件 检测 引物 对 
的 退火 温度 Tm CG 含量 ,以 及 是 否 会 产生 二 聚 体 ,利用 OligoCale 软件 检测 引物 自身 的 互补 性 ,包括 潜在 的 发 
夹 结构 .3' 端 的 互补 性 、 自 身 退 火 位 点 的 个 数 等 基本 参数 。 
1.2 引物 敏感 性 、 特 异性 评估 

采用 在 线 软件 TaxMan ,根据 核糖 体 基因 数据 库 , 以 多 来 源 的 rDNA 为 模板 ,通过 多 序列 全 局 比 对 ,对 引物 
扩 增 的 生物 类 和 群 和 未 扩 增 的 生物 类 群 物种 数量 进行 预测 统计 及 比较 ,进而 对 3 对 引物 的 敏感 性 NS 特异 性 进行 
评估 。 本 实验 选取 已 发 表 的 海洋 真 核 生物 高 通 量 测序 引物 V9(F/AR) 和 C4(F/R) 作 为 比较 分 析 参 照 物 "9 。 
13 浮游 植物 鉴定 
1.3.1 形态 学 鉴定 

分 别 在 辽东 湾 长 兴 岛 海域 采集 3 份 浮游 植物 水 样 ,经 福 尔 马 林 固定 后 带 回 实验 室 进 行 显微镜 检测 。 
13.2 ”分子 鉴定 

(1) 样 品 采集 

现场 采集 1L 表层 海水 ,经 20 um 微 孔 滤 膜 过 滤 去 除 小 型 及 大 型 浮游 生物 ,然后 用 0.22 pm 微 孔 滤 膜 过 滤 
收集 微型 和 微微 型 浮游 生物 ,最 后 将 滤 膜 转移 至 1.5 mL 无 菌 离心 管 中 , 钾 于 -20 % 或 -80 CARET .运输 。 

(2)DNA 提取 

采用 CTAB 法 提取 浮游 生物 基因 组 ,将 0.22 um 滤 膜 剪 碎 于 了 5 mL 离心 管 中 , 加 入 500 pL CTAB 裂解 液 
(2% CTAB;100 mmol/L Tris-Cl pH 8.0;1.4 mmol/L NaCl; 10 mmol/L EDTA) I 1 uL B- 琉 基 乙 醇 ,5 一 10 uL 4& 
FE K,559C 裂解 1— 1.5 h。 短 暂 离心 ,取出 液体 于 新 的 离心 管 中 , 用 酚 握 仿 抽 提 2 次 后 , 取 上 清 ,加 入 两 倍 体 
积 预 冷 的 无 水 乙醇 ,沉淀 2 一 3 h ,保留 沉淀 ,使 用 ,75% 乙 醇 清洗 沉淀 ,得 到 微型 浮游 生物 基因 组 DNA,1% 琼 脂 
糖 凝 胺 电泳 检测 并 利用 紫外 分 光 光 度 计 测定 WBNA 浓度 及 纯度 ,-20% 冰箱 保存 备用 。 

(3) 引 物 合成 和 PCR 扩 增 

引物 序列 设计 完成 , 交 由 上 海 生 开 生物 公司 进行 合成 。PCR 反应 体系 为 50 kL, 包括 PCR Buffer 5 pL、 
dNTP Mixture 8uL 、 上 下 游 引 物 (10 kmolL) 各 2 kL 模板 DNA 2 pL,Taq DNA 聚合 酶 2.5 U, 加 适量 灭 菌 水 。 
扩 增 反应 均 在 PE 9700 型 了 CR [E CPE 公司 ) 上 完成 ,反应 条 件 :94% 预 变性 3 min;94% 变性 30 s,58% 退火 45 
s,72*C 延伸 45 s, 483 个 循环 ;72% 延伸 5 min。1% 琼 脂 糖 凝 胶 电 泳 检 测 PCR 产物 。 将 检测 合格 产物 交 由 诺 
和 致 源 生 物 信 息 科 技 有 限 公 司 , 使 用 NEB Next@ Ultra™ DNA Library Prep Kit for Ilumina ( New England 
Biolabs) 建 库 试 剂 盒 进 行文 库 的 构建 ,构建 好 的 文库 经 过 Qubit 定量 和 文库 检测 ,合格 后 ,使 用 Hiseq2500 
PE250 HERE HLIEF - 

(人 数据 分 析 

测序 得 到 的 原始 数据 使 用 FLASH 软件 进行 拼接 ,参照 Qiime 软件 质量 控制 流程 ,将 拼接 后 的 序列 经 过 截 
取 、 过 滤 得 到 有 效 数据 。 利 用 Uparse 软件 对 有 效 数据 进 OTUs( Operational Taxonomic Units) 聚 类 和 物种 分 类 
分 析 ,RDP Classifier 方法 与 Silva 数据 库 对 OUTs 代表 序列 进行 物种 注释 分 析 。 再 对 OTUs 进行 丰 度 .多样 性 
此 数 等 分 析 ,同时 对 物种 注释 在 各 个 分 类 水 平 上 进行 群落 结构 的 统计 分 析 。 最 后 在 以 上 分 析 的 基础 上 ,可 以 
进行 一 系列 的 基于 OTUs 和 物种 组 成 的 聚 类 分 析 ,挖掘 样品 内 和 样品 间 的 物种 组 成 差异 。 


2 结果 


2.1 引物 设计 
利用 软件 MEGAA 对 179 条 序列 进行 多 重 序列 比 对 ,找到 对 应 的 V4 和 V9 区 ,同时 结合 测序 技术 要 求 进 
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行 引 物 设计 ,片段 长 度 为 300 一 450bp, 并 使 用 FPCR ,OligoCale 等 软件 对 引物 的 基本 参数 进行 检测 和 评估 ,其 
基本 参数 满足 引物 设计 的 基本 要 求 ,引物 设计 结果 见 表 1。 


表 1 微型 藻类 18S rDNA 可 变 区 扩 增 引物 序列 


Table 1 Primers used in the 18S rDNA variable region of microalgae 


引物 名 称 Primer 引物 序列 Sequence( 5'-3') 引物 名 称 Primer 引物 序列 Sequence( 5'-3') 

V4-F 5'-GCGGTAATTCCAGCTCCAATA-3' V9-R 5'-CCTTCYGCAGGTTCACCTAC-3' 
V4-R 5'-GATCCCCHWACTTTCGTTCTTGA - 3' C4-F 5'-CCAGCASCYGCGGTAATTCC-3' 
V9-F 5'-CCCTGCCHTTTGTACACAC-3' C4-R 5'-ACTTTCGTTCTTGATYRA-3' 


214 引物 敏感 性 特异 性 评估 
采用 在 线 软 件 Taxman 对 3 对 引物 的 敏感 性 .特异 性 进行 评估 ,结果 如 表 2。 统 计数 据 最 示 s 引 物 V9 (F 
R) 对 的 整体 扩 增 效率 低 于 引物 V4(F/R) 和 CA C EAR) ,引物 V4(F/R) 对 真菌 .后 生动 物 \ 领 懂 毛 虫 目 类 群 的 


扩 增 特异 性 低 于 CACE/R) ,对 定 甘 藻 纲 \ 绿 藻 门 等 生物 类 群 的 扩 增 特异 性 高 于 CACF/R) ,ULB] VA(F/R) 在 扩 
增 真 核 藻类 类 群 中 较 有 优势 。 


表 2 引物 对 扩 增 SILVA 数据 库 中 真 核 生 物 类 群 的 特异 性 统计 
Table 2 Specificity of the PCR primers foreukaryotarDNA sequences in the SILVA database 


分 类 单元 真 核 生物 序列 eukaryotic sequence 分 类 单元 真 核 生物 序列 eukaryotic sequence 

Taxon V9( F/R) V4(F/R) CA(F/R) | Taxon V9( F/R) V4( F/R) C4( F/R) 
真菌 Fungi 0.11(816) 0.29(2165) 0.42(3122) HEI ]/Dinophyceae 0.19(178) — 0.61(583)  0.68(650) 
后 生动 物 Metazoa 0.12(1524) ”0.44(5628) 0.70(8900) || Katel J] Cryptophyta 0.28(37) 0.91 (121) 0.88(117) 
SiR H Choanoflagellida ”0.36(16) 0.68(30) 0.89(39) | FERETE Dictyochophyceae 0.40(8) 0.95(19) 1.00(20) 
植物 界 Viridiplantae 0.17(691)  0.75(3137) 0.75(3115) | 43 Chrysophyceae 0.29(42) | 0.86 (123)  0.85(122) 
AE SIEGE] Haptophyceae 0.32(46) | 0.80(114) — 0.01 (2) || &ESET] Bacillariophyta 0.23(150) 0.77 (514) — 0.82(547) 
ktel] Chlorophyta 0.27(344) | 0.70(888) | 0:67(850) | 黄 群 藻 纲 Synurophyceae 0.28(20) 0.72 (52) 0.72(52) 


引物 对 主要 扩 增 真 核 类 群 的 频数 分 布 ,( ) 内 的 数字 表示 扩 增 该 类 群 的 序列 数 ,( ) 前 的 数字 表示 扩 增 序列 数 占 该 类 群 序列 数 的 比例 


2.2 引物 对 基因 组 DNA 扩 增 产物 检测 

为 了 比较 3 对 引物 在 实际 样品 电 对 真 核 藻类 扩 增 特 
异性 的 差异 ,对 采集 的 3 份 样品 分 别 进行 了 引物 扩 增 。 扩 
增产 物 的 检测 结果 见 图 1“ 本 文中 3 了 3 对 引物 VA(F/R) 、V9 500bp 


500bp 
(F/R) FI C4( F/R) KID 385779] Fr BEA 30] 2g 370bp ,135bp, 、 250pp 250bp 
= S yo x 100b 
380bp 左右 ,图 1-9, Fr Be ERR UI i R R E E 250bý 
符 , 扩 增产 物 量 满足 测序 需求 ,PCR 体系 合格 。 
2.3 不同 引 物 对 的 测序 结 VA(F/R) V9(F/R) CA(F/R) 
38 T: HluminaHiSeq2500 测序 平台 PE250 测序 ,得 到 的 1 引物 扩 增 片段 检测 


原始 数据 使 用 FLASH 软件 进行 拼接 ,参照 Qiime 软件 质量 Fig.1 PCR amplification of representative DNA templates 
控制 流程 ,将 拼接 后 的 序列 经 过 截取 、 过 滤 得 到 有 效 数 据 。 with rDNA primer sets 
每 人 样品 平均 获得 68834 条 原始 序列 ,经 过 拼接 和 质量 过 
滤 , 每 个 样品 平均 得 到 68420 条 有 效 序列 ,高 质量 数据 占 到 99% 以 上 ,获得 基因 注释 的 序列 数 在 97% 以 上 , 表 
明 测 得 的 数据 准确 可 靠 。9 个 样品 测序 结果 见 表 3。 

为 了 研究 样品 的 物种 组 成 多 样 性 信息 ,使 用 Uparse 软件 对 所 有 样品 的 有 效 序列 进行 聚 类 ,以 97 多 的 一 臻 
性 将 序列 聚 类 成 为 OTUs。 每 个 样品 平均 获得 308 个 OTUs ,每 个 OTUs 代表 一 类 物种 。3 对 引物 V4( F/R), 
V9(F/R) 和 C4(F/R) 鉴 定 的 平均 微型 /微微 型 浮游 植物 OTUs 数 分 别 为 78、42 .58。 
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表 3 有效 序列 数 及 OTUs 数 
Table 3 The statistics of Effective Tags and OTUs numbers 


引物 名 称 样品 名 称 原始 序列 有 效 序列 注释 序列 总 OUT 种 水 平平 均 OTUs “种 水 平 浮游 植物 OTUs 
Primer Sample Name Raw Tags Effective Tags Taxon Tags Total OUT num Average OUT num Phytoplankton OUT num 
V4( F/R) LDW1.1 61435 61213 57690 254 145 68 
LDW2.1 65094 64871 61954 272 148 73 
LDW3.1 57105 56319 55193 323 176 92 
V9( F/R) LDWI.2 63565 63547 59444 293 62 37 
LDW2.2 62681 62674 60858 461 74 46 
LDW3.2 77199 77177 75234 434 69 44 
C4(F/R) LDWI.3 64345 64117 63330 151 98 33 
LDW2.3 73653 73198 71258 250 139 59 
LDW3.3 93829 92663 90861 334 191 81 


2.4 镜 检 优 势 种 属 与 分 子 鉴 定 结果 比较 


镜 检 优 势 种 属 


与 高 通 量 测序 结果 见 表 4。 从 结果 中 可 以 看 出 引物 V4( F/R) 对 镜 检 优势 种 属 的 检 出 率 与 
引物 C4( F/R) 相 当 , 并 且 都 高 于 引物 V9( F/R) ,对 于 样品 中 物种 多 样 性 的 挖掘 ,引物 W4(F/R) 略 高 于 引物 


C4(F/R), 且 明显 高 于 引物 V9(F/R)。 


mA 镜 检 优势 种 属 与 分 子 检测 结果 


Table 4 The microscopic examination.results of three samples 


3 NT m 镜 检 DES MOST 
"uf 样品 编号 种 类 Microscopic examination The numberofalgae species 
z Sample Species names poem 
密度 Density/ ( M/L) VA(F/R) V9(F/R) C4( F/R) 
: = ) 1 Thalassiosira sp. 3.6 4 1 3 
j Chaetoceros sp. 183.6 2 1 1 
2 Chaetoceros sp. 606.7 1 4 1 
Noctilucascintillans 563.3 1 0 1 
3 Chaetoceros sp. 102.5 2 0 2 
- Asterionellopsisglacialis 315.0 1 0 1 
) 根据 表 5 注释 结果 ,在 种 水 平 下 ,统计 3 对 引物 对 应 的 测序 鉴定 结果 ,引物 VACE/R) 鉴定 出 的 浮游 植物 


数量 多 于 V9 (F/R 


AF C4 (E/R) 鉴定 出 的 数量 。3 对 引物 对 优势 种 Dolichomastix tenuilepis, 细小 微 胞 藻 
( Micromonas pusilla DN 4x "E flt ER SE ( Ostreococcus tauri ) 、 um pk a ( Pycnococcus provasolii ) 、 里 氏 金 色 藻 
( Chrysochromu linaleadbeateri ) AE AERE ( Aureococcus anophagefferens ) JR SES E ( Heterosigma akashiwo ) 等 


微型 和 微微 型 浮游 植物 的 鉴定 能 力 比较 ,引物 VACE/R) 与 C4(F/R) 相 近 , 二 者 强 于 VO(F/R) 。 


表 5 3 对 引物 鉴定 浮游 植物 种 类 
Table 5 Amplification of Phytoplankton with three rDNA primer pairs 


物种 名称,Species names 拉丁 文 Latin names V4(F/R) V9( F/R) C4( F/R) 
HE] Phylum Dinoflagellata Akashiwo sp. + 
Alexandrium pseudogoniaulax + 
Alexandrium tamarense + 
Amoebophrya sp. + + 
Amylax triacantha $ + 
Ceratium tenue + 
Duboscquella sp. + + 
Gonyaulax polygramma + + 
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续 表 

物种 名 称 Species names 拉丁 文 Latin names V4(F/R) V9( F/R) C4(F/R) 
Gyrodinium sp. 4 " 
Karlodinium sp. 十 
Karlodinium veneficum 十 
Neoceratium sp. 十 
Noctiluca scintillans + + $ 
Pelagodinium beii + 
Prorocentrum micans + + $ 
Prorocentrum minimum i 
Protoperidinium bipes + 
Protoperidinium pellucidum + 
Scrippsiella sp. + 4 
Scrippsiella trochoidea E 
Stoeckeria. sp. " 
Stoeckeria. algicida 十 $ 
Woloszynskia sp. + 4 
Woloszynskia halophila + 

绿 藻 门 Phylum Chlorophyta Bathycoccus sp. 十 d 
Crustomastix sp. T E 十 
Dolichomastix tenui 十 + 
Mamiella gilva + + + 
Mantoniella squamata E 
Micromonas pusilla 十 十 4 
Micromonas sp. 十 
Nannochloris sp. 十 十 十 
Nephroselmis pyriformis 十 E E 
Ostreococcus tauri + + + 
Picochlorum sp. 
Prasinophyceae sp. RCC39 * 十 + 
Pseudococcomyxa simplex 十 
Pterosperma cristatum 十 + 
Pyenococcus provasolii + + E 
Pyramimonas sp. 十 十 
Tetraselmis sp. 十 E 

patel] Phylum Cryptophyceae Geminigera cryophila + 
Hemiselmis sp. + 4 
Hemiselmis cryptochromatica + + 
Katablepharid sp. 十 + 
Leucocryptos sp. + E 
Leucocryptos marina 十 十 E 
Myrionecta rubra 十 
Rhodomonas sp. E 
Teleaulax gracilis 十 + 
Teleaulax sp. + E 
Teleaulax amphioxeia 十 十 d 

XE WESEL] Phylum Haptophyta Calyptrosphaera sp. LKM-2007-1 十 十 
Chrysochromulina sp. 十 + 
Chrysochromu linaleadbeateri + + 
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续 表 
物种 名 称 Species names 拉丁 文 Latin names V4(F/R) V9( F/R) C4( F/R) 
Emiliania sp. + + 
Haptolina fragaria + + 
Imantonia sp. + 
Prymnesiumn palpebrale 十 
Tergestiella adriatica + 
EET Asterionellopsis glacialis + + + 
Phylum Bacillariophytina Chaetoceros calcitrans + + + 
Chaetoceros muellerii $ + + 
Chaetoceros socialis + 
Chaetoceros sp. p442 + F + 
Cyclotella choctawhatcheeana + 
Ditylum brightwellii + + + 
Entomoneis cf. alata + + 
Hyalosira sp. t 
Minutocellus sp. $ + 
Peridiniumquin quecorne endosymbiont + 
Pleurosigma planktonicum + 
Rhizosolenia setigera + 
Rhizosolenia similoides + 
Skeletonema sp. * 十 
Skeletonema pseudocostatum 十 
Stephanopyxis turris 十 十 
Thalassiosira sp. 十 十 
Thalassiosira concaviuscul, 十 十 
Thalassiosira mala 十 
Thalassiosira minima 十 十 
Thalassiosira-oceanica F $ k 
KEME] Apedinella radians 十 十 十 
Phylum Heterokontophyta Aureocoecus anophagefferens 十 十 十 
Bolidomonas pacifica 十 
Dictyocha fibula 十 
Ochromonas sp. t 十 
Dictyocha speculum + + + 
Florenciella parvula + + + 
Heterosigma akashiwo + F F 
Paraphysomonas sp. + + 
Pedinella sp. squamata + 
Phaeocystis antarctica + 
Pseudochattonell afarcimen 十 十 
Pseudopedin ellaelastica 十 十 十 
Pteridomonas sp. t 十 
合计 64 57 57 


T 


+ 表示 引物 扩 增 片段 的 注释 结果 在 3 个 样品 中 任意 出 现 一 次 即 标记 + 


3 讨论 


分 子 生物 学 技术 极 大 地 推动 了 微型 真 核 营 类 多 样 性 的 研究 。 目 前 ,世界 各 大 海域 的 微微 型 浮游 生物 分 子 
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多 样 性 逐步 开展 研究 ,相关 基因 文库 数据 也 在 不 断 更 新 和 扩充 ?1。 针 对 海洋 微型 和 微微 型 藻类 鉴定 方面 
存在 的 难题 ,本 实验 选择 18S rDNA 的 可 变 区 作为 目标 基因 ,构建 高 通 量 测序 文库 ,应 用 生物 信息 学 处 理 数据 ， 
对 海洋 微型 和 微微 型 藻类 进行 了 多 样 性 研究 。18S rDNA 结构 和 功能 保守 ,区 分 力 可 达 种 属 水 平 , 常 用 于 真 核 
生物 物种 鉴定 和 系统 发 生 学 分 析 。18SrDNA 全 长 序列 交替 排列 着 保守 区 和 可 变 区 ,反映 了 生物 种 间 的 亲缘 关 
系 和 物种 间 的 差异 ,利用 了 保守 区 设计 引物 , 扩 增 出 可 变 区 ,根据 不 同 种 属 可 变 区 的 序列 特征 进行 分 类 :2 。 
真 核 生物 共 含 有 8 个 可 变 区 ,其 中 V4 和 V9 区 是 目前 广泛 使 用 的 真 核 生物 标记 基因 中 。Wuyts 等 比较 分 析 
了 3253 条 核糖 体 rDNA 核酸 序列 的 V4 区 ,发 现 该 区 域 为 真 核 生 物 核糖 体 rDNA RI AERE PECES Dod 
Stoeck 和 Micah 等 人 研究 发 现 ,V4 区 在 区 分 亲缘 关系 较 近 的 物种 具有 较 大 优势 :9 , 目 与 V9 区 相 比 ,V4 区 得 
到 的 实验 结果 更 接近 于 以 18S rDNA 全 长 为 目标 基因 的 研究 结果 2 。18S rDNA 的 V4 区 变 蜡 较 大 长度 较 
长 .均一 性 好 ,适用 于 高 通 量 测 序 技术 ,是 微微 型 藻类 鉴定 的 良好 目标 基因 。 

本 文通 过 对 3 对 引物 的 敏感 性 特异 性 进行 评估 结果 显示 :引物 VAER 对 真菌 Ji ^E 2147 AE d H 
类 群 的 扩 增 特异 性 要 低 于 V9( E/R) fI CA(E/R) ,而 对 真 核 藻类 的 扩 增 特异 性 均 高 于 9GEBMR ) , 扩 增 红 藻 门 
的 特异 性 低 于 CACE/R) 。 为 了 验证 引物 VA(F/R) , V9( FR) 和 C4( F/R) 对 海水 梢 品 中 微型 /微微 型 浮游 杆 
物 多 样 性 分 析 的 可 行 性 和 差异 性 ,采用 TlluminaHiSeq2500 测序 平台 PE250 测序 , fet both EPIKE FP EXSI] 
物 对 应 的 平均 微型 /微微 型 浮游 植物 OTUs 数 , 并 与 样品 的 显微镜 检测 结果 进 条 比较 zj 结果 显示 ,引物 对 V4 
(F/R) 在 微型 /微微 型 浮游 植物 种 类 数 鉴定 方面 优 于 其 他 两 对 引物 。 高 通 量 测 疗 获得 浮游 植物 种 类 与 显微镜 
观察 结果 进行 比较 ,引物 V4( F/R) 对 镜 检 结果 中 的 藻类 的 检 出 率 与 引物 C4(F/R) 相 同 ,并 且 都 高 于 引物 V9 
(F/R) ,对 于 样品 中 物种 多 样 性 的 挖掘 能 力 , 引 物 VA (E/R) ESTE 519] C4(F/R), 且 明显 高 于 引物 V9(F/ 
R) 。 

根据 注释 结果 ,在 种 水 平 下 ,统计 3 对 引物 对 应 的 扩 增 产物 的 鉴定 结果 ,3 对 引物 对 优势 种 细小 微 胞 藻 、 
AE ORBE ARE ERRER . 抑 食 金 球 菠 .未 漳 异 变 薄 等 微型 和 微微 型 浮游 植物 的 鉴定 能 力 ,引物 V4(F/ 
R) 与 CA(F/R) flit, — T. V9(F/R) 。 其 中 ,金牛 微 球 藻 是 一 种 单 细胞 绿 洛 ,直径 0.8um 左右 ,是 目前 发 
现 的 个 体 最 小 .细胞 结构 简单 基因 组 结构 却 十 分 复杂 的 自 养 真 核 细胞 生物 '”] Vaquer 557 AUT 3S ESL PE 18 
ill ( Thaulagoon ) 金牛 微 球 藻 丰 度 与 铀 含量 有 关 , 并 发 现 贝 类 无 法 摄食 这 种 藻类 。 赤 潮 异 弯 藻 是 一 种 细胞 微 
小 .分布 广 泛 . 广 温 、 广 盐 性 赤潮 藻 种 ,对 楚 他 生物 影响 较 大 '”。 抑 食 金 球 灌 , 呈 金光 色 球 形 或 亚 球 形 ,直径 
2pm 左右 ,属于 微微 型 藻类 ,是 目前 发 现 引 发 褐 潮 的 主要 种 类 ,对 生态 环境 和 汐 业 养殖 影响 较 大 "3。2013 年 
7 月 在 长 兴 岛 近 岸 海域 出 现 术 人 异常 现象 ,经 分 子 检 测 , 主 要 优势 种 为 抑 食 金 球 藻 , 丰 度 超过 10' 个 /L' 引 。 于 
杰 等 人 利用 通用 引物 扩 增 18S vDNA. 可 变 区 V9 后 ,结合 高 通 量 技术 ,建立 浮游 生物 多 样 性 高 效 检 测 技术 ,对 
渤海 海域 的 4 个 褐 潮 多 发 区 的 优势 藻 种 进行 分 析 , 发 现 抑 食 金 球 藻 和 多 形 微 眼 藻 共同 引发 了 该 海域 的 
18197, 

本 研究 表明 LV/ 188 rDNA 的 V4 区 作为 目标 基因 ,设计 引物 V4( F/R) ,借助 IlluminaHiSeq 2500 PE250 高 
TEFAF AMERRE ,建立 了 海洋 微型 和 微微 型 藻类 高 效 检测 方法 。 该 方法 既 传 水 了 传统 分 子 技 
术 的 特点 汉 具 入-“ 些 独特 的 优势 :测序 通 量 高 ,实验 流程 简单 ;数据 处 理 快速 .准确 ;灵敏 度 高 ;节约 成 本 ; 适 
用 性 广泛 针对 海洋 中 不 同 水 域 .水 层 时间 的 微型 /微微 型 浮游 植物 多 样 性 的 研究 均 适 用 。 


参考 文献 ( References) : 


[1] BEWE. 厦门 港 微型 浮游 硅 藻 研究 [D]. 厦门 : 厦门 大 学 , 1990. 

[2] LeeSR, Oak J H, Chung I K, Jin A L. Effective molecular examination of eukaryotic plankton species diversity in environmental seawater using 
environmental PCR, PCR-RFLP, and sequencing. Journal of Applied Phycology, 2010, 22(6) : 699-707. 

[3] Balzano S, Marie D, Gourvil P, Vaulot D. Composition of the summer photosynthetic pico and nanoplankton communities in the Beaufort Sea 
assessed by T-RFLP and sequences of the 18S rRNA gene from flow cytometry sorted samples. The ISME Journal, 2012, 6(8) : 1480- 1498. 
[4] WuL, Yu Y, Zhang T L, Feng W S, Zhang X, Li W. PCR-DGGE Fingerprinting Analysis of Plankton Communities and Its Relationship to Lake 
Trophic Status. International Review of Hydrobiology, 2009, 94(5) : 528-541. 


http ://www.ecologica.cn 


201706.00849v1 


" 
LI 


chinaXiv 


ChinaXiv 合 作 期 刊 


4216 E dm 学 R 3735 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 
[20] 


[21] 
[22] 
[23] 
[24] 
[25] 


[26] 


[27] 


Handelsman J, Rondon M R, Brady S F, Clardy J, Goodman R M. Molecular biological access to the chemistry of unknown soil microbes; a new 
frontier for natural products. Chemistry &Biology, 1998, 5(10) : R245-R249. 

Massana R, Karniol B, Pommier T, Bodaker I, Béjà O. Metagenomic retrieval of a ribosomal DNA repeat array from an uncultured marine 
alveolate. Environmental Microbiology, 2008, 10(5) : 1335-1343. 
Simon C, Daniel R. Metagenomic analyses; past and future trends. Applied and Environmental Microbiology, 2010, 77(4) : 1153-1161. 

陈 随 清 , 王 利 丽 . rRNA 基因 (rDNA) 序 列 分 析 在 中 药品 种 鉴定 中 的 应 用 及 研究 进展 . 河南 中 医药 学 报 ，2003，8(2) : 86-88. 

Neefs J M, Van de Peer Y, Hendriks L, De WachterR. Compilation of small ribosomal subunit RNA sequences. Nucleic Acids Research, 1990, A8 
(Suppl) : 2237-2317. 


Van Der Staay S YM, De Wachter R, Vaulot D. Oceanic 18S rDNA sequences from picoplankton reveal unsuspected eukaryotic diversity. Nature, 
2001, 409(6820) : 607-610. 

López-García P, Rodríguez-Valera F, Pedrós-Alió C, Moreira D. Unexpected diversity of small eukaryotes in deep-sea Antarctic plankton. Nature, 
2001, 409(6820) : 603-607. 

BS, xii. KHZ, JERES, WIEK. 自然 水 样 微型 藻类 遗传 多 样 性 的 方法 学 研究 . 海洋 科学 ,2002，26(4) "1-4. 

xi, RBS, KAF, RAE, JERES. 南沙 群岛 微型 与 超 微型 真 核 藻 类 遗传 多 样 性 的 初步 研究 . 海洋 科学 ,2003, 27(7) : 43-47. 

Yuan J, Chen M Y, Shao P, Zhou H, Chen Y Q, Qu L H. Genetic diversity of small eukaryotes from the coastal waters of Nansha Islands in China. 
FEMS Microbiology Letters, 2004, 240(2) : 163-170. 

Sogin M L, Morrison H G, Huber J A, Welch D M, Huse S M, Neal P R, Arrieta J M, Herndl G J. Microbial diversity in the deep sea and the 


[zz 


underexplored " rare biosphere” . Proceedings of the National Academy of Sciencesof the United States of America, 2006, 103(32) : 12115-12120. 
Stoeck T, Bass D, Nebel M, Christen R, Jones M D, Breiner H W, Richards T A. Multiple marker parallel tag environmental DNA sequencing 
reveals a highly complex eukaryotic community in marine anoxic water. Molecular Ecology, 2010, 19(81) : 21-31. 

Behnke A, Engel M, Christen R, Nebel M, Klein RR, Stoeck T. Depicting more accurate pictures of protistan community complexity using 
pyrosequencing of hypervariable SSU rRNA gene regions. Environmental Microbiology, 2011, 13(2) : 340-349. 

Amaral-Zettler LA, McCliment E A, Ducklow H W, Huse S M. A method for studying protistan diversity using massively parallel sequencing of V9 
hypervariable regions of small-subunit ribosomal RNA genes. PLoS One, 2009, 4(7) : e6372. 

Worden A Z. Picoeukaryote diversity in coastal waters of the Pacific Ocean. Aquatic Microbial Ecology, 2006, 43(2) : 165-175. 

Marie D, Zhu F, Balagué V, Ras J, Vaulot D. Eukaryotic picoplankton. communities of the mediterranean sea in summer assessed by molecular 
approaches (DGGE, TTGE, QPCR). FEMS Microbiology Ecology, 2006, 55(3) : 403-415. 

李 绿 砚 . 中 国 南海 微型 与 超 微型 真 核 浮 游 生物 的 分 也 生态 党 研究 [DD]. 广州 : 中 山大 学 , 2008. 

江 雪 娇 . 北 黄海 微微 型 浮游 植物 的 丰 度 及 微微 型 起 核 浮 游 生物 分 子 多样 性 研究 [D]. 青岛 : 中 国 海洋 大 学 ,2009. 

Pawlowski J, Christen R, Lecroq B, Bachar D, Shahbazkia HR, Amaral-Zettler L, Guillou L. Eukaryotic richness in the abyss: insights from 
pyrotag sequencing. PLoS One, 2011, 6(4) : e18169. 


KI J S. Hypervariable regions ( V1-V9 yof the-dinoflagellate 18S rRNA using alarge dataset for marker considerations. Journal of Applied Phycology , 
2012, 24(5) : 1035-1043. 

Wuyts J, De Rijk P, Van de,Peer Y,-Pison G, Rousseeuw P, De Wachter R. Comparative analysis of more than 3000 sequences reveals the 
existence of two pseudoknots in area. V4 of eukaryotic small subunit ribosomal RNA. Nucleic Acids Research, 2000, 28(23) : 4698-4708. 
Dunthorn M, Klier J} Bunge J, Stoeck T. Comparing the Hyper-Variable V4 and V9 regions of the small subunit rDNA for assessment of ciliate 
environmental diversity. Journal of Eukaryotic Microbiology, 2012, 59(2) : 185-187. 

Courties C, Perasso R, Chrétiennot-Dinet M J, Gouy M, Guillou L, Troussellier M. Phylogenetic analysis and genome size of Ostreococcus tauri 
( Chloropliyta, Prasinophyceae). Journal of Phycology, 1998, 34(5) : 844- 849. 

VaquerA, TroussellierM, CourtiesC, BibentB. Standing stock and dynamics of picophytoplankton in the Thau Lagoon ( northwest Mediterranean 
coast)» Limnology and Oceanography, 1996, 41(8) : 1821-1828. 
Bü, ATL, 钱 培 元 . 赤潮 异 弯 藻 Heterosigma akashiwo 的 生长 特性 . 海洋 与 湖沼 ，2002,，33(2) :209-214. 

Bowers H A, Tomas C, Tengs T, Kempton J W, Lewitus A J, Oldach D W. Raphidophyceae [ Chadefaud Ex Silva V] systematics and rapid 
identification; sequence analyses and Real-Time PCR Assays. Journal of Phycology, 2006, 42(6) : 1333-1348. 


Keppler C J, Hoguet J, Smith K, Ringwood A H, Lewitus A J. Sublethal effects of the toxic alga Heterosigma akashiwo on the southeastern oyster 
( Crassostrea virginica). Harmful Algae, 2005, 4(2) : 275-285.. 

Probyn T, Pitcher G, Pienaar R, Nuzzi R. Brown tides and mariculture in Saldanha Bay, South Africa. Marine Pollution Bulletin, 2001, 42(5) : 
405-408. 

辽宁 省 海洋 与 渔业 厅 . 2013 年 辽宁 省 海洋 环境 状况 公报 . 中 国 海洋 报 , 2014-04-09. 

于 杰 . 浮游 生物 多 样 性 高 效 检测 技术 的 建立 及 其 在 渤海 褐 淹 研 究 中 的 应 用 [D]. 青岛 : 中 国 海洋 大 学 , 2014. 


http ://www.ecologica.cn 


